Urban Traffic Dynamics: A Scale-Free Network Perspective 
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This letter propose a new model for characterizing trafiic dynamics in scale-free networks. With 
a replotted road map of cities with roads mapped to vertices and intersections to edges, and intro- 
ducing the road capacity L and its handling ability at intersections C, the model can be applied to 
urban traffic system. Simulations give the overall capacity of the traffic system which is quantified 
by a phase transition from free flow to congestion. Moreover, we report the fundamental diagram 
of flow against density, in which hysteresis is found, indicating that the system is bistable in a 
certain range of vehicle density. In addition, the fundamental diagram is significantly different from 
single-lane traffic model and 2-D BML model with four states: free flow, saturated flow, bistable 
and jammed. 
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PACS numbers: 89.40.-a, 45.70.Vn, 89.75.Hc, 05.70.Fh 



Traffic and transportation is nowadays one of the 
most important ingredients of modern society. We rely 
greatly on networks such as communication, transporta- 
tion and power systems. Ensuring free traffic flow on 
these networks is therefore of g reat research interest 
[lllllSlllllillllllSllJ Recently, more and 
more empirical evidence indicates that these networked 
systems are of small-world and scale-free structural fea- 
tures 

Cmillllll And so the traffic flow on scale-free 
networks is being widely investigated 0, 0, 0, 0, ■ 
In the present work, we propose a new model for the traf- 
flc dynamics of such networks. The potential application 
of this research will be urban (road) traffic networks. 

Previously work on urban traffic normally maps roads 
to edges and intersections to vertices. In 1992, Biham, 
Middleton and Levine (BML) 8] proposed an ideal 2- 
dimensional model for studying urban traffic. They used 
a, N X N grid to represent the urban transportation net- 
work. Initially, cars are placed at the sites. Each car is 
either East-facing or North-facing. At odd time steps, 
each North-facing car moves one unit North if there is 
a vacant site before it. At even time steps. East-facing 
cars move East in the same way. The rules of BML model 
can be considered as controlling urban traffic by traffic 
lights. The BML model reproduced the transition from 
free-flow phase to jammed phase. Since then, many stud- 
ies have been carried out on the base of the BML model. 
For example, Nagatani P investigated the effect of car 
accident induced jam on the phase transition of urban 
traffic; Chung et al. ^3 investigated the effect of im- 
properly placed traffic lights; and recently. Angel et al. 
[ll| discussed the jammed phase of the BML model. 

The models mensioned above always adopt a, N x N 
grid (or close to that) to represent the urban traffic sys- 
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FIG. 1: Illustration of the road network mapping. Each 
straight road in (a) is mapped to a vertex in (b). And inter- 
sections are mapped to edges between the vertices. 



tem. However, real urban traffic system is obviously 
much more complicated. Perhaps, the most natural way 
is to map each intersection to a vertex and each segment 
of road to an edge between vertices. However, obviously, 
this kind of simulation is computation-consuming and the 
interaction of any two neighboring segments of a straight 
road is in most cases wiped off. 

The unique feature of our model is that we look at the 
urban traffic networks in a different point of view. We 
will create traffic flow networks with roads mapped to 
vertices and intersections to edges between vertices, as 
was inspired by the information networks suggested in 
12] for an information view of a city. In this way(See 
Fig^, the degree k of vertex is the number of intersec- 
tions along the street and a major road with many minor 
roads connected to it can be seen as a major vertex. Em- 
pirical observations demostrate that the remapped urban 
networks exhibit scale free structural properties. 

With this new paradigm, we can look at urban traffic 
networks from a novel perspective. In network language, 
a trajectory of a vehicle can then be mapped in an urban 
traffic network map from a road (vertex) to another road 
(vertex) through a one directional channel (edge) of an 
intersection. In this work we will take this perspective 
to the extreme, and assume that the travel tiine/cost of 
just driving along a given road can be zero jl^l- Our 
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model is partially inspired by the work of information 
packet flow on the Internet 111 El El El El . The 
phase transition from free flow to jammed state, and the 
plot of flow against density that all have been empirically 
observed can be well reproduced by this model. 

Though the previous work found that the replotted 
road networks of cities have scale- free characteristics E2 j 
there is no well-accepted model for these networks up to 
now. Without losing generality, our simulation starts 
from generating the urban transportation network ac- 
cording to the most general Barabasi- Albert scale-free 
network model Q|. In this model, starting from mo fully 
connected vertices, one vertex with m edges is attached 
at each time step in such a way that the probability Hi of 
being connected to the existing vertex i is proportional to 
the degree ki of the vertex, i.e. 11^ = '^here j runs 

over all existing vertices. The capacity of each vertex 
(road) is controlled by two parameters: (1) its maximum 
cars number L, which is proportional to its degree k (a 
long road ordinarily has more intersections and can hold 
more cars): L = a x k; (2) the maximum number of 
cars handled per time step, which reflects the capability 
of intersections: C = /3 x L. Motivated by the Internet 
information flow models 11 El El El Ea , the system 
evolves in parallel according to the following rules: 

1. Add Cars - Cars are added with a given rate R (cars 
per time step) at randomly selected vertices and each car 
is given a random destination. 

2. Navigate Cars - If a car's destination is found in 
its next-nearest neighborhood, its direction will be set to 
the destination vertex. Otherwise, its direction will be set 

to a neighboring vertex h with probability: Pfi = . 

Here the sum runs over the neighboring vertices, and (p is 
an adjustable parameter. It is assumed that the cars are 
unaware of the entire network topology and only know 
the neighboring vertices' degree fc^. 

3. Cars Motion - At each step, only at most C cars 
can leave a vertex (road) for other vertices and FIFO 
(first-in-first-out) queuing discipline is applied at each 
vertex. When the queue at a selected vertex is full, the 
vertex won't accept any more vehicles and the vehicle 
will wait for the next opportunity. Once a car arrives at 
its destination, it will be removed from the system. 

We first simulate the traffic on a network of = 100 
vertices (roads) with mO = m = 2, a = 5 and f3 = 0.2. 
This relatively small system can be seen as simulating the 
backbone of a city's urban traffic network. The selection 
of (3 is based on the single-road traffic flow theory which 
shows that the maximum flux on a highway is about 20% 
of its maximum density [3, H H ■ For simplicity, we 
do not consider the phenomenon that the flux decreases 
when the density is above 20%. The combination of a 
and P can be also interpreted as: each intersection can 
handle one car turning for one road at each step. 

To characterize the system's overall capacity, we first 
investigate the car number Nc increment rate u> in the 
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FIG. 2: (color online). The overall capacity of a road network 
with iV = 100, mO = m = 2, a = 5 and /3 = 0.2. (a) The 
variation of car number Nc for different R when (p = 0.1. 
Rc{= 13) is determined at the point where the Nc increment 
rate ui increases suddenly from zero and Nc increases rapidly 
towards the system's maximum car number, (b) The critical 
Rc versus 4>. The maximum of Rc corresponds to (/> = 0.1 
marked by a dash line. The data are obtained by averaging 
Rc over 10 network realizations. 



At) — Nc{t)) takes average over time windows of width 
At. Fig|2Ia) shows the variation of Nc with different 
R for (j) — 0.1. One can see that there is a critical Rc 
(= 13) at which Nc runs quickly towards the system's 
maximum car number and uj{R) increases suddenly from 
zero. uj{R) = corresponds to the cases of free flow state, 
which is attributed to the balance between the number 
of added and removed cars at the same time. However, 
if R exceeds the critical value Rc, cars will in succession 
accumulate in the system and then congestion emerges 
and diffuses to everywhere. Ultimately almost no cars 
can arrive at their destinations. 

Evidently, Rc is the onset of phase transition from free 
flow to jammed state. Hence, the system's overall ca- 
pacity can be measured by the critical value of Rc under 
which the system can maintain its normal and efficient 
functioning. Fig|2fb) depicts the variation of Rc versus 
(j). The maximum overall capacity occurs at = 0.1 
(shghtly greater than 0.0) with i?™"^ = 13. 

Here we give a heuristic analysis for determining the 
optimal value of (jj with the maximum capacity. If we 
neglect the queue length L of each vertex, for = 0, cars 
will move in the system nearly like random walk. There 
is a well-known result from graph theory that if a particle 
performs a random walk, in the limit of long times, the 
time the particle spends at a given vertex is proportional 
to the degree k of the vertex EH ■ One can easily find out 
that, the number of cars observed averagely at a vertex 
is proportional to the degree of that vertex. Meanwhile, 
the cars handling ability of each vertex is assumed to be 
proportional to its degree. Thus in the case oi (j> — 0, this 
rule produces an average effect that no congestion occurs 
earlier on some vertices with particular degree than on 
others. Accordingly, (j) — results in the maximum sys- 
tem capacity. However, in our model, each vertex has a 
limited queue length L — a x k and R cars are gener- 
ated randomly among all vertices at each time step, so 
small degree vertices are slightly more easily congested. 
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FIG. 3: (color online). Average travel time (T) versus </> for 
R = 1 and 2. The data are truncated because the system 
jams when (j) is either too large or too small. The right panel 
shows the variation of (T) versus R when is fixed. The data 
are also truncated when the system jams. 



Therefore, for our traffic model, a (p slightly larger than 
zero can enhance the system's capacity maximally. 

Then we simulate the cars' travel time spent in the ur- 
ban transportation system. It is also an important fac- 
tor for measuring the system's efficiency. In FigO^a) , we 
show the average travel time (T) versus (j) under traffic 
load B. = 1 and 2. In the free-flow state, almost no con- 
gestion on vertices occurs and the time for cars waiting in 
the vertex queue is negligible, therefore, the cars' travel 
time is approximately equal to their actual path length 
in replotted road map. But when the system is close 
to a jammed state, the travel time will increase rapidly. 
One can see that when (j) is close to zero, the travel time 
is minimum. In Fig|SJb) inset, the average travel time is 
much longer when (f) is negative than it is positive. These 
results are consistent with the above analysis that a max- 
imum Rc occurs when (j) is slightly greater than zero. Or, 
in other words, this effect can also be explained as fol- 
lows: when > 0, cars are more likely to move to the 
vertices with greater degree (main roads), which enables 
the main roads to be efficiently used and enhance the sys- 
tem's overall capability; but when (f> is too large, the main 
roads will more probably get jammed, and the efficiency 
of the system will decrease. 

Finally, we try to reproduce the fundamental diagram 
(flux-density relation) of urban traffic system. It is one 
of the most important criteria that evaluates the transit 
capacity for a traffic system. Our model reproduced the 
phase transition and hysteresis in fundamental diagram. 

To simulate a conservative system (constant density), 
we count the number of arriving cars at each time step 
and add the same number of cars to randomly selected 
vertices of the system at the beginning of next step. The 
flux is calculated as the number of successful car turn- 
ings from vertex to vertex through edges per step, as is 
similar to the Internet information flow. Here we ignore 
the movement of cars on a given road. In fact, the flux of 
car turnings at intersections can, to some extent, reflect 
the flux on roads. In Fig^ the fundamental diagrams 
for (j) ~ 0.1,0.0, 1.0 and —0.5 are shown. The curves of 
each diagram show four flow states: free flow, saturate 
flow, bistable and jammed. For simplicity, we focus on 
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FIG. 4: (color online). Fundamental diagram for & N = 100 
network with mO — m — 2, a = 5, /3 = 0.2, and different 
(j). The data are averaged over 10 typical simulations on one 
realization of network. In each chart, the solid square line 
shows the flux variation when adding cars randomly to the 
system (increase density), while the empty circle line shows 
the flux variation when drawing out cars randomly from the 
system (decrease density). And the data are collected and 
averaged at 10,000-11,0000 steps after the intervention when 
the system has reached a steady state. The sudden transition 
density values are: 0.76 and 0.45 ((/) = 0.1), 0.82 and 0.76((?!) = 
0.0), 0.26 and 0.22 {(p = 1.0), 0.83 and 0.80 ((^ = -0.5). For 
different realiazations of network, the charts are similar in 
phases, but with different transition values. 



the = 0.1 chart in the following description. As we can 
see, when the density is low (less than w 0.1), all cars 
move freely and the flux increases linearly with car den- 
sity. It is the free-flow state that all vertices (roads) are 
operated below its maximum handling ability C. Then 
the flux's increment slows down and the flux gradually 
comes to saturation (0.10 ~ 0.45). In this region, the flux 
is restricted mainly by handling ability C of vertices. One 
can see that when (j) is close to zero, the saturated flux 
(~ 360) is much higher than other values. 

At higher density, the model reproduces an important 
character of traffic flow - "hysteresis". It can be seen 
that two branches of the fundamental diagram coexist 
between 0.45 and 0.76. The upper branch is calculated 
by adding cars to the system, while the lower branch is 
calculated by removing cars from a jammed state and 
allowing the system to relax after the intervention. In 
this way a hysteresis loop can be traced (arrows in Fig^. 
The hysteresis loop indicates that the system is bistable 
in a certain range of vehicle density. And as we know 
so far, it is the flrst time that our model reproduces the 
hysteresis phenomenon in scale-free network traffic and 
in urban network traffic. 

To test the finite-size effect of our model, we simu- 
late some bigger systems with much more vertices(roads). 
The simulation shows similar phase transition and hys- 
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FIG. 5: (color online), (a) Fundamental diagram for a, N — 
1000 network with mO = m = 5, a = l,/3 = 0.2 and 4> = 0.1. 
(b) The averaged number of jammed vertices (Njv)- The 
symbols for increasing/decreasing density are the same as in 
FigH One can see that the two sudden change points 0.32 
and 0.21 in both charts are equal. 

teresis in fundamental diagram as shown in Fig[3Ja). 

The flux's sudden drop to a jammed state from a 
saturated flow is a first order phase transition. This 
behavior can be explained by the sudden increment of 
full(jammed) vertices in the system (See FiglSJb)). Ac- 
cording to the evolution rules, when a vertex is full of 
cars, the cars at neighboring vertices can not turn to 
it. So the cars may also accumulate on the neighboring 
vertices and get jammed. This mechanism can trigger 
an avalanche across the system when the car density is 
high. As shown in Fig[Sl the number of full vertices in- 
crease suddenly at the same density where the flux drop 
to zero and almost no car can reach its destination. As for 
the lower branch of the bistable state, starting from an 
initial jammed configuration, the system will have some 
jammed vertices that are difficult to dissipate. Clearly, 
these vertices will decrease the system efficiency by af- 
fecting the surrounding vertices until all vertices are not 
jammed, thus we get the lower branch of the loop. 



Moreover, an important conclusion can be drawn by 
comparing the = 0.1 chart with the (j) = 0.0 chart in 
FigElthat the — 0.1 chart has a much broader bistable 
region than the — 0.0 one. This means, when the 
system retreats from a heavy load jammed state, it is 
more difficult to reach a high efficiency state if 4> is greater 
than zero that cars are more likely to move to main roads. 
In other words, though it is wise to take full advantage 
of the main roads when the entire traffic is light, it won't 
be so happy to do so at rush hours. 

In conclusion, a new traffic model for scale-free net- 
works is proposed. In the new perspective of mapping 
roads to vertices and intersections to edges, and incorpo- 
rating road/intersection capability limits, the model can 
be applied to urban traffic system. In a systemic view 
of overall efficiency, the model reproduces several signifi- 
cant characteristics of network traffic, such as phase tran- 
sition, travel time, and fundamental diagram. A special 
phenomenon - the "hysteresis" - can also be reproduced. 
Although the microscopic dynamics of urban traffic are 
not well captured, such as the movement of cars along 
streets, the interaction with traffic lights, and the dif- 
ferences between long and short streets, our model is 
still a simple and good representation for urban traffic, 
since much empirical evidence is well reproduced by the 
model. Further effort is deserved to consider more de- 
tailed elements for better mimicking real traffic systems. 
Moreover, by choosing other values of the parameters, 
the model may be applied to other networked systems, 
such as communication and power systems. 
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